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What is the order of the deconfining phase transition? 

A. Di Giacomo a 

a Dipartimento di Fisica, Universita di Pisa, 
I.N. F.N. Sezione di Pisa 

The techniques are discussed by which the order of the deconfining phase transition is investigated on the 
lattice. QCD with Nf = 2 is a special case, which can provide information on the mechanism of confinement. 
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1. Introduction 

Quarks and gluons have never been observed 
as free particles, a phenomenon known as color 
confinement. 

A widely accepted idea in the 1960's was 
that a limiting temperature exists for hadronic 
matter 1 , the Hagedorn temperature Tjj; in- 
creasing the energy density just below Th would 
result in production of particles, rather than in 
an increase of temperature. In ref.|2] it was sug- 
gested that a phase transition could exist to a 
plasma of quarks and gluons (Deconfining Tran- 
sition) . 

Big experiments to explore this possibility in 
heavy ions collisions have not yet given conclusive 
results 0. 

The only existing source of information are 
Montecarlo simulations of the theory on the Lat- 
tice. 

The partition function of a system of fields at 
temperature T 



Z = Tr{exp(-H/T)} 



(1) 



is equal to the Euclidean Feynmann integral of 
the theory extending in time from to 1/T, with 
periodic boundary conditions in time for bosons, 
antipcriodic for fcrmions. In the Lattice formu- 
lation this corresponds to simulate the theory 
on a lattice x N t , with the spatial extension 
N s 3> Nt, the time extension. The temperature 
is T = l/a(0,m q )Nt, where (3 = 2N/g 2 and a is 
the lattice spacing in physical units. 

To explore the deconfining transition a signa- 
ture for confincmcnt-dcconfincmcnt must be de- 



fined. In quenched QCD (Nf — , no dynami- 
cal quarks) this is done as follows. Consider the 
Polyakov line 



L{x) — Pexp 



The correlator 




A (x,t)dt 



(2) 



G(r) = (L(f)£t(0)) (3) 
is related to the static qq potential by the relation 



V{r) = -T\nG{r) 



(4) 



On the other hand at large distances cluster prop- 
erty requires 



G(f) 



cexp(-^r) + |(L)| 2 



T 



(5) 



When (L) = cq© gives a linearily rising po- 
tential at large distances, V = err which is taken 
as a definition of confinement. When (L) ^ 0, 
V(r) — > const and there is no confinement. (L) 
is an order parameter for confinement, the corre- 
sponding symmetry being Z3. 

A deconfining phase transition is indeed ob- 
served in lattice simulations, corresponding to the 
above definition , both in SU(2) and in SU(3) 
pure gauge theories: (L) grows rapidly from to 
1 by increasing T at some T c . In a finite lattice , 
however, the growth of (L) is smooth: no phase 
transition can take place in a finite system 0]. 
The raise of (L) becomes sharper and sharper as 
the volume increases. 

The infinite volume limit can be studied by a 
technique used in statistical mechanics [S], known 
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as finite size scaling analysis. The susceptibility 
of the order parameter xl can be defined at a 
given spatial size N s 



Xl(N s 



cPx(L(x)tf(0) - L(0)i f (0)) (6) 



Xl gives a measure of the slope of the increase of 
(L(x)) at T c . As N s — + oo it diverges with some 
critical index 7 when approaching T c from below 



T 

a- E ) 



(7) 



The correlation length £ of the order parameter 
(L(x)) also diverges at the critical point with a 
critical index ^ 

t * r- v (8) 

T— 

Assuming that the ratio 

is an analytic function, as T — > T c , a/£ — > and 
the scaling law follows 



i?(iV s )c/(0,^) 



or 



(9) 



(10) 



The dependence on A s is dictated by the critical 
indices of the transition. In particular the value 
of xl at the peak r = scales as 

XL,peak(N s ) CX N2 /U (11) 

A similar discussion gives for the specific heat C v 

C v = C° v + N?' v *c{NI' v t) (12) 

with a the corresponding critical index and C® 
an additive term due to the presence of an ad- 
ditive renormalization. By studying numerically 
the dependence of xl and C v on N s the critical 
indices a, 7 can be determined and the order 
of the transition and the universality class with 
them. A weak first order transition is a limiting 
case when a = 7 = 1 and v = 1/d with d the 
number of spatial dimensions. 



For quenched SU(2) it is found that the uni- 
versality class is that of the 3d ising model, and 
v = .62jHj. For quenched SU(3) the transition is 
weak first order and v = 1/3. 

An alternative order parameter (/x) is the vev 
of a magnetically charged operator [S]. This oper- 
ator will be described in some detail in the next 
section for full QCD. A finite size scaling analysis 
of 

p = -^H») 

gives results consistent with those obtained by 
using (L) and (fj,) ^ for T < T c , (/x) = for 
T > T c . 

Since (//) 7^ signals dual superconductivity, 
this result can be considered as evidence that dual 
superconductivity of the vacuum is the mecha- 
nism of confinement. 

2. Full QCD 

In the presence of dynamical quarks a clear way 
to define and to detect confinement does not exist. 
(L) is not a good order parameter, since the cor- 
responding symmetry Z3 is broken by the quark 
coupling. In addition string breaking is expected 
to occur: the static potential energy converts into 
dynamical q — q pairs at large distances so that 
the potential is not linear any more with the dis- 
tance, even if there can be confinement. 

At m q = chiral symmetry exists, which is 
spontaneously broken at T = 0, the pseudoscalar 
mesons being the corresponding Goldstone par- 
ticles and the quark condensate (tptp) the or- 
der parameter. At some critical temperature 
T c ~ 170 Mev chiral symmetry is restored. It 
is not well understood what is the interplay of 
chiral symmetry with confinement. Moreover at 
m q 7^ chiral symmetry is explicitely broken and 
(V'V') i s n °t a good order parameter either. 

However if one looks at susceptibilities like 
X£i (7« they all show in numerical simula- 
tions a maximum for a T c (m q ) as functions of T, 
along a line in the plane (m q , T) , which is as- 
sumed by convention as a critical line , with the 
confined phase below it and the deconfincd phase 
above itPfTU]. 
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Figure 1. Phase diagram of two flavor QCD 



We shall discuss a simplified model with Nf — 
2 and m u = rrid = to, which is semirealistic , but 
is of special interest as a probe of the theory. The 
critical line is schetched in fig 1. for this model. 

A renormalization group analysis of the chiral 
transition at m = leads to the following pre- 
dictions , based on the hypothesis that the Gol- 
stone particles are the relevant degrees of freedom 
at the transition. For Nf = 3 or 2 + 1 the tran- 
sition is first order, and such is expected to be 
for For N* = 2 the transition is first or- 

der both at to = and at to if the axial U(l) 
broken by the anomaly is restored at a lower tem- 
perature than chiral symmetry,it is second order 
if the anomaly persists up to the chiral transition, 
and in that case the line in fig. 1 is a crossover . 
There exists no precise numerical analysis of the 
order of the transition , but for some reason the 
second possibility is usually accepted. 

One can assume dual superconductivity of the 
vacuum as a criterion for confinement. This cri- 
terion works in the quenched case, where it is 
consistent with the one based on (L) [H]. To be 
rigorous one should prove that dual superconduc- 
tivity implies confinement of any colored particle 
. In the same way, however, one should prove 
that (L) = implies confinement of any colored 
particle. 

Dual superconductivity is detected by a disor- 
der parameter {/J, a ), the vev of an operator carry- 



ing non zero magnetic charge[Hj-A non zero vev of 
[i a signals Higgs breaking of magnetic U(l) sym- 
metry, or dual superconductivity. We shall recall 
below the definition and some properties of \x a . 
Numerical simulations show indeed that vacuum 
is a dual superconductor in the region of fig 1 be- 
low the critical line , and is normal above it|H] 

The operator /i a is defined 12 by the orbits 
U(x)cf>d iag U^(x) in the gauge group of the sym- 
metric space defined by the co-roots 



diag 



= di 




The definition is 

U(x) 



(14) 
abelian 



Each choice of U(x) defines an 
proiection|13|. and with it monopoles. The oper- 
ator fi a creates a monopole in the given abelian 
projection. It can be shown[T3 | |T5 | . however, that 
the property of (/i a ) of being non-zero or zero is 
abelian projection independent, if the number 
density of monopoles is finite. Numerical simu- 
lations show that this is indeed the case. If one 
looks at the distribution of the difference of eigen- 
values of operators in the adjoint representation 
on the lattice sites, the number of sites in which 
that difference vanishes is zero. That difference 
is equal to zero at the locations of monopoles. A 
large sample of configurations, lattice spacings, 
and operators have been studied. 

A finite size scaling analysis allows to deter- 
mine the critical indices and the order of the 
transition^] The problem has two scales, the 
correlation length and the quark mass so that 
neglecting the ratio of the lattice spacing to the 
correlation length, 



<M>=/(f, 



iN2») 



(15) 



where yh is the anomalous dimension of to . By 
choosing values of to and N s such that the second 
argument is constant a scaling law follows for the 
susceptibility 
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whence v can be extracted. The result is consis- 
tent with a first order thansition [y = 1/3) 17 . 

This determination should be consistent with 
the scaling of the specific heat eq(12). Prelimi- 
nary data show that this is the case indeed [TT]. 
Further numerical work is on the way to put the 
result on a firm basis. This would definitely legiti- 
mate the assumption of (/x) as an order parameter 
for confinement, and confirm dual superconduc- 
tivity as a mechanism for color confinement. It 
would also imply 11 that cither axial U(l) sym- 
metry is restored before chiral symmetry, or that 
Goldstone particles are not the relevant degrees of 
freedom at the critical point. The first possibility 
can also be tested on the lattice. 
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